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ABSTRACT 



We study the signatures of coronal heating on the differential emission mea- 
sure (DEM) by means of hydrodynamic simulations capable of resolving the 
chromospheric-corona transition region sections of multi-stranded coronal loops 
and following their evolution. We consider heating either uniformly distributed 
along the loop or localized close to the chromospheric footpoints, in both steady 
and impulsive regimes. Our simulations show that condensation at the top of 
the loop forms when the impulsive heating, with a pulse cadence lower than 
the plasma cooling time, is localized at the loop footpoints and the pulse en- 
ergy is below a threshold above which the heating balances the radiative losses, 
thus preventing the catastrophic cooling which triggers the condensation. A con- 
densation does not produce observable signatures in the DEM because it does 
not redistribute the plasma over a sufficiently large temperature range. On the 
other hand, the DEM coronal peak is found sensitive to the pulse cadence time 
when this is longer or comparable to the plasma cooling time. In this case, 
the heating pulses produce large oscillations in temperature in the bulk of the 
coronal plasma, which effectively smears out the coronal DEM structure. The 
pronounced DEM peak observed in active regions would indicate a predominance 
of conditions in which the cadence time is shorter or of the order of the plasma 
cooling time, whilst the structure of the quiet Sun DEM suggests a cadence time 
longer than the plasma cooling time. Our simulations give an explanation of 
the warm overdense and hot underdense loops observed by TRACE, SOHO and 
Yohkoh. However, they are unable to reproduce both the transition region and 
the coronal DEM structure with a unique set of parameters, which outlines the 
need for a more realistic description of the transition region. 
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Introduction 



Solving the problem of the heating of the solar corona is one of the major issues in solar 
physics. In particular, the question whether the plasma heating inside coronal structures is 
the effect of steady or impulsive processes, uniform or localized within the structure, is still 
open. 

Since the earliest observations, the solar corona appeared as composed by discrete 
bright structures, called coronal loops, consisting of magnetic flux tubes in which the hot 
and dense plasma is confined by the magnetic field. Coronal loops are characterized by 
different lengths, temperatures, activity levels, and appear t o evolve with lifetimes of the 
order of several hours. The first models of loop heating (e.g. 



Rosner et al 



1978 



Serio et al. 



198ll ) considered flux tubes of constant cross-section filled with plasma in hydrostatic 
equilibrium, and in energy balance under the effects of uniform steady heating, conductive 
heat flux, and radiative losses. These models predict scaling laws relating the temperature, 
density, and length of a loop, and allow us to reproduce quite satisfactorily the X-rays 
coronal emission of both the Sun and solar-type stars. 

Nevertheless, recent TRACE and SOHO observations have provided evidence that a 
large majority of warm coronal lo ops (T ~ 1-2 MK) , alth o ugh appearing in quasi-static 



conditions, are indeed o ver dense (lAschwanden et al. 



Patsourakos et al. 



1999 



2001 



Winebarger et al. 



2003 



20041 ). while hot loops (T > 2 MK) observed by Yohkoh exhibit lower 



densi ties than predicted by loop models in hydrostatic equilibrium (IPorter fc Klimchuk 



19951 ) . These discrepancies could be explained if coronal loops are assumed to consist 
of unresolved magnetic strands, ea ch of them heated impulsively and non-uniformly, at 



different times from its neighbors ( Car gill 



2003 



2006 



Cargill fc Klimchuk 



Klimchuk et al 



2004 



Reale et al. 



994 



2005 



Klimchuk fc Car gill 



2001 



Patsourakos fc Klimchuk 



Spadaro et al. 



2005 



Klimchuk 



20081 ). The idea of impulsive coronal heating was firstly proposed by 
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Parkerl (119831 . Il988l ). who introduced the concept of nanoflare, a local small-scale (< 100 km) 
event in which roughly 10 24 erg of energy are released in the corona by magnetic field 
reconnection, after the magnetic stresses induced on the flux tubes by photospheric motions 
have reached a certain threshold. 

In this Paper we investigate the response of magnetic loop plasma to different kinds 
of energy deposition by performing numerical simulations with a code capable of resolving 
the chromospheric-coronal transition region sections of the loop and following them as they 
respond to heating variations. We consider heating either uniformly distributed along the 
loop or localized close to the chromospheric footpoints, in both steady-state and impulsive 
regimes. In the impulsive case, we consider different values for the cadence of the injection 
of the energy pulses into the coronal segment of the loop, in order to perform a first 
comprehensive study of the consequences of the variation of such a parameter for the loop 
structure and evolution. The hydrodynamic behavior of coronal loops undergoing different 
heating regimes and the relevant consequences on the differential emission measure are 
discussed in detail, together with the indications that the variety of conditions found in this 
exploration give on the physical origins of coronal heating and related phenomena. 



Numerical model and simulations 



We performed numerical simu 



heating conditions using ARGOS (lAntiochos et al. 



ations of corona' 



loop 



1999 



hydrodynamics under different 



MacNeice et al. 



20001 ) . a one- 



dimensional code that solves the standard set of equations for the conservation of mass, 
momentum and energy by means of a high-order Godunov scheme and an adaptive mesh 
refinement. ARGOS also properly accounts for the heat fl ux saturation th at may occur in 



some low-density and high-temperature regimes (see, e.g. 



Klimchu 



20061 ) 
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We considered a flux tube with a 80 Mm long coronal segment, to model typical 



active region loops, as observed, for instance, by TRAC 
characteristic features of the loop model can be found in 



A detailed descri ption of the 



Spadaro et al. 



fl2003h . 



The loop was initially settled to a nearly stable equilibrium state under the action 
of a spatially uniform and t emporally constant background heating (as described, e.g, in 



Antiochos et al 



1999 



2000j). A heating rate -Ebase = 2 x 10 5 erg cm 3 s 1 was chosen 



according to the canonical hydrostatic scaling laws (see 



Rosner et al. 



19781 ). to get an apex 



temperature of about 0.75 MK at the e nd of the relaxation phase. At this stage, the loop 



structure is that of a steady loop as in 



Rosner et al. 



(119781 ). The uniform background 



heating applied to achieve the initial equilibrium was turned off, within the coronal segment 
of the loop, when we started the additional heating (t = 0) and remained off for the 
remaining of the simulation. In the chromospheric part of the loop, the E^ ase heating 
was maintained throughout the simulation, although it has a little influence on the loop 
hydrodynamics, due to the high density in the chromosphere. 

The additional heating rate, E(s,t), is assumed to be a separable function of the 
curvilinear coordinate along the field lines, s, and of the time, t: 

E(s,t) = F(s) ■ G(t) [erg cm" 3 s" 1 ] . (1) 

We examined all the possible combinations among two kinds of heating spatial distribution, 
i.e. quasi-uniforno or localized at the loop footpoints, and two kinds of temporal 
distribution, i.e. impulsive or steady. 

In the localized heating cases, the heating rate has a maximum at the loop footpoint, 
and exponentially falls off in the corona with a fixed scale-length of 10 Mm; the location 



x Note that in the quasi-uniform case the chromospheric and coronal heating rates are 
different. 
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and scale-len gth of the energy deposit i on are consistent with those deduced from TRACE 



observations (lAschwanden et al. 



2000 



F(s) 



20011 ). We also consider the possibility of an 
asymmetric energy deposition between the two footpoints of the loop, since real coronal 
loops do not appear to be symmetric. For the impulsive heating cases, we inject into 
the coronal segment of the loop a sequence of energy pulses, or nanoflares, with constant 
cadence, duration, and energy amplitude; we model each nanoflare with a Maxwellian 
function, such that the energy release E(s,t) during a single event has a steep rise, followed 
by a more gradual decrease. Specifically: 

/exp (— ^x 1 ) , s > Si localized heating 
1 quasi-uniform heating 

(2) 

. £/2^t 2 exp (-£) impulsive heating 

Es steady heating 

where / = 0.75 or 1.0 measures the constant ratio of the localized heating at the right 
footpoint to that at the left one, s\ is the initial position of the top of the chromosphere 
at each footpoint, A = 10 Mm is the scale-length of the heating distribution, Ej gives the 
volumetric heating per nanoflare event in the impulsive case, r is a characteristic parameter 
related to the duration of the energy pulses (here always equal to 12.5 s), and E s is the 
constant heating rate in the steady case. Table [1] lists the values of the parameters adopted 
in the simulations performed for this study. 

We chose the value of Ej so that 

f dV f E(s,t)dt = E P [erg], (3) 
Jv Jo 

where V is the volume of the coronal section of the loop, t c is the cadence time between 
two consecutive nanoflares, and Ep is the amount of energy supplied by each pulse. For the 
steady heating cases, the value of E$ was adjusted so that the total energy deposited into 
the loop in a time interval t c was equal to Ep, as listed in Tabled! In this case Es is equal 
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Table 1: Parameters of the simulations 



Run 


A 




Hip 


E b 


<E> c 


Class 


# 


(Mm) 


(s) 


(10 24 erg) 


(erg cm -3 s _1 ) 


(10 21 erg s" 1 ) 




1 


Quasi-uniform 


Steady 


1 


5 x 10" 4 


4 


Steady 


2 


Quasi-uniform 


250 


1 


2 x 10~ 3 


4 


Steady 


3 


Quasi-uniform 


500 


1 


2 x 10~ 3 


2 


Dynamic 


4 


Quasi-uniform 


1000 


1 


2 x 10~ 3 


1 


Dynamic 


5 


Quasi-uniform 


2000 


1 


2 x 10~ 3 


0.5 


Dynamic 


6 


10 


Steady 


1 


4 x 10~ 3 


4 


Condensation 


7 


10 


250 


0.125 


2.75 x 10~ 3 


0.5 


Condensation 


8 


10 


250 


0.5 


1.1 x 10~ 2 


2 


Condensation 


9 


10 


250 


1 


2.2 x 10~ 2 


4 


Condensation 


10 


10 


250 


2 


4.4 x 10~ 2 


8 


Steady 


11 


10 


250 


4 


8.8 x KT 2 


16 


Steady 


12 


10 


500 


0.25 


5.5 x KT 3 


0.5 


Dynamic 


13 


10 


500 


1 


2.2 x 10~ 2 


2 


Dynamic 


14 


10 


500 


2 


4.4 x 10~ 2 


4 


Dynamic 


15 


10 


1000 


1 


2.2 x 10~ 2 


1 


Dynamic 


16 


10 


1000 


4 


8.8 x 10~ 2 


4 


Dynamic 


17 


10 


2000 


1 


2.2 x 10~ 2 


0.5 


Dynamic 


18 


10 


2000 


8 


1.76 x 10- 1 


4 


Dynamic 



a Total energy per pulse, or, in the steady cases, total energy deposited into the loop in a time interval equal 
to t c = 250 s. 

b Maximum heating rate per unit volume. 

c Average amount of energy supplied to the loop per unit time. 
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to E max , as reported in Tabled! 



The value of the time interval between the pulses, i.e. the cadence t c , was fixed taking 
into account the characteristic radiative cooling time of the loop, r coo i- For a loop with a 
semilength L = 40 Mm and an apex temperature T apex = 2 MK, it results r coo i 1000 s, 
according to, e.g., 



Serio et al. 



(1199 lh . Here we report results obtained with t c ranging from 



250 s « r cool /4 to 2000 s « 2r cool . 

We also considered different values for Ep, in order to change the amount of energy 
supplied to the loop by the sequence of heating pulses and investigate the related effects on 
the plasma hydrodynamics. 



3. Results and discussion 

3.1. Plasma dynamics 

Figure [1] shows the initial part of the temporal evolution of the plasma temperature, 
density, a nd velocity averaged over a re 



(see, e.g. 



Patsourakos fc Klimchuk 



. evant portion (3/4) of the loop coronal segment 
20051 1 in some representative cases. 

In the impulsive quasi-uniform heating case with Ep = 10 24 erg and t c = 250 s (run 2), 
the loop model simply settles into a new quasi-static equilibrium state, characterized by 
nearly constant temperature, density and velocity, apart from small oscillations, lower then 
10%, due to the sequence of energy pulses. Note that the average values of these parameters 
behave as the corresponding ones obtained for the steady quasi-uniform heating case 
(run 1). After the start of the additional heating, the temperature of the loop immediately 
reaches a value of about 2.3 MK; then, in the next hour, it decreases and stabilizes around a 
value slightly lower than 2 MK. The density monotonically increases from 1.4 x 10 8 cm~ 3 to 
~ 9.0 x 10 s cm~ 3 , while the plasma velocity, after a first highly dynamic phase caused by the 
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Fig. 1. — Initial part of the temporal evolution of the model temperature (top panels), 
density (middle panels), and velocity (bottom panels), averaged over the upper 3/4 of the 
loop coronal segment, for run 1, 2, and 5 (left panels; black, heavy gray and dotted lines, 
respectively), and run 6, 9, 10, and 17 (right panels; black, heavy gray, light gray, and dotted 
lines, respectively). The velocity is almost identical in runs 9 and 10 and therefore only the 
run 9 velocity is reported in the bottom-right panel. The vertical dashed lines in the right 
panels approximatively correspond to the maximum of the condensation phase for run 9, 
whose snapshot is given in Fig. [2j 
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abrupt increase of the heating rate and characterized by upflow velocities up to 40 km s l , 
approaches very low values of some km s _1 . The time-averaged plasma temperature and 



density for t 



lis loop model are consistent with the hydrostatic scaling laws defined by 



Rosner et al. 



(119781 ) . 



Conversely, asymmetric heating localized near the loop footpoints, both steady (run 6) 
and impulsive with t c = 250 s (run 9), causes a more dynamic evolution for similar values 
of the heating parameters, with long-term fluctuations of temperature and density (after 
an initial steep increase) due to cycles of plasma condensation formation, motion along the 
loop, and falling onto the nearest and less heated footpoint (note that the first cycle only is 
reported in Fig. [p. 



The phenomenon of plasma condensation formation is well known anc 



described in the literature (e.g., 



Karpen fc Antiochos 



2008 



Antiochos et al. 



Klimchuk et al 



2009 



1999 



2000 



parpen et a 



Miiller et al 



2003 



2004 



extensively 



2001 



2006 



Testa et al 



O'Shea et al. 



2005). Recently, it has also been observed in coronal non-flaring loops (e.g. 
20071 ). It is the effect of a thermal instability occurring near the top of the loop, where the 
energy supply is not sufficient to balance the radiative losses induced by the increase of 
the plasma density owing to the strong chromospheric evaporation. This eventually leads 
to a catastrophic cooling, down to chromospheric temperatures, of the material located 
close to the loop apex and, consequently, to the formation of a region of low-temperature 
and high-density plasma. This is clearly visible in Fig. [2], which shows the instantaneous 
plasma temperature and density profiles along the loop (run 9) at a stage of its evolution 
characterized by the presence of a condensation. The blob of cool and dense plasma 
subsequently starts moving slowly toward the less heated footpoint, because of the 
differences in pressure induced by the asymmetric energy release in the two legs of the loop, 
and finally drains onto the chromosphere. The depleted loop then reheats quickly and a new 
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2.5 F 




10 8 " 



40 60 80 100 120 140 160 

s (Mm) 



Fig. 2. — Instantaneous profiles of the loop model temperature (upper panel) and density 
(lower panel) vs. the curvilinear coordinates for run 9, at 2.7 hours after the start of the 
additional heating, marked with a vertical line in Fig. [TJ The chromospheric segments have 
been truncated to highlight the coronal portion of the loop. 
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cycle of chromospheric evaporation, plasma condensation and draining takes place. In the 
present model the first catastrophic cooling phase approximatively begins 2.5 hours after 
the start of the additional heating and lasts for about one hour; its total duration slightly 
depends on whether the energy deposition is impulsive or steady, in the latter case being 
nearly 20 minutes shorter. The periodicity of the described cycle is about three hours. 

The temporal evolution of run 7 and 8 (with Ep < 10 24 erg, not shown in Fig. [1]) is 
very similar to that of run 9, with, however, lower temperatures, densities, and velocities, 
but again with the onset of a dynamic cycle of plasma condensation formation. 

Therefore, the impulsive heating cases, both quasi-uniform (run 2) and localized 
(run 7, 8, and 9), do not exhibit appreciable differences with respect to the corresponding 
steady cases (run 1 and 6, respectively). This is expected when the cadence of nanofl ares 



is very rapid compared to t 



Testa et al. 



2005 



Klimchuk 



re lo op plasma cooling time (see, e.g. 



Walsh et al. 



1997 



20061 ). In our cases a cadence time equal to r coo i/4 already 



gives a nearly steady heating situation. On the other hand, it appears evident that the 
localization of the heating near the loop footpoints plays a fundamental role in the plasma 
condensation formation when the energy deposition is steady or impulsive with a cadence 
time well below the loop cooling time. These results, obtained for a ratio of the heating 
damping-length (A) to the loop semilength equal to 1/4, are in agreem ent with those of 



other works modeling 



20041 : 



Testa et al. 



2005|). 



ootpoint-heated loops in similar conditions (e.g. 



Miiller et al 



2003 



Run 5 and 17 illustrate both quasi-uniform and localized heating cases in which 
the cadence time is longer than the loop plasma cooling time (t c > r coo i). In this cases 
(dotted-line curves in Fig. [I]) the temperature in the upper part of the loop coronal segment 
shows pronounced oscillations since there is enough time between pulses for the plasma to 
cool down and drain downward to the chromospheric region (as confirmed by the downward 
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velocities of the plasma noticed in between the pulses). Density has small variations, but 
velocity shows strongly damped oscillations with a rather high peak at the beginning of 
each energy pulse. No condensation forms in this case, because the intermittent heating 
prevents the accumulation of plasma at the loop top, and thus the thermal instability. A 
similar behavior is also obtained when t c ~ r coo i (run 4, 15, and 16) or equal to r coo \/2 
(run 3, 12, 13, and 14). 

Note, however, that although the heating localization and a high frequency pulse 
cadence are necessary to yield a catastrophic cooling phase during the loop evolution, they 
are not always sufficient. The crucial point is the balance between the energy supplied 
to the loop top by the sequence of heating pulses and the radiative losses of the plasma 
accumulated therein. In fact, by sufficiently increasing the amount of energy supplied by 
each pulse, even without changing the other parameters of the heating regime, such a 
balance can be achieved, thus preventing the thermal instability and the consequent plasma 
condensation formation. This is the case of run 10 (see left panels in Fig. [1]) and run 11. 

As far as the quasi-uniform loop simulations (runs 1-5) are concerned, it is also worth 
noting that increasing t c between runs 2-5 corresponds to a decrease by a factor of 8 in the 
total average energy <E> put into the corona per second, because E P is held fixed. Since 
no condensation forms, the average temperature and density are essentially controlled by 
the mean energy dissipated per unit time, as in the steady case (see also Sect. 13.2.2]) . 

Figure [3] reports the average loop density and temperature during the evolution of 
some representative loops together with the hydrostatic scaling law relationship. We note 
that all the examined loops exhibit densities significantly smaller than those predicted 
by the scaling laws in the initial phase of their evolution (for about one hour), in both 
steady-state and impulsive cases. Even if the onset of the heating produces an increase in 
the loop temperature and a large conductive heat flux at the coronal base of the loops, 
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Fig. 3. — Evolution of the average loop density and temperature for runs 1, 2, and 5 (top 
panel; solid, dash-dotted, and dotted line, respectively), run 6, 9, and 17 (middle panel; 
solid, dash-dotted, and dotted line, respectively), and run 10 and 11 (bottom panel; dash- 
dotted and dashed line, respectively). The hydrostatic scaling law (heavy dashed line, see 



Rosner et al. 



19781 ) is also reported, gauged according to the initial static equilibrium (in- 



dicated by the cross). In the middle panel, the evolution characterized by condensation 
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driving the upward expansion known as chromospheric evaporation (lAntiochos fc Sturrock 



19781 ). there is a time lag in the corresponding increase of the coronal density, so that the 



loops appear to be hot and underdense. 

The subsequent evolution of the loop depends on the characteristics of the plasma 
heating. 

In the quasi-uniform cases (top panel in Fig. [3]), when the heating is steady (run 1) or 
impulsive with a cadence time well below the loop cooling time (run 2), the loop evolves 
toward a new state almost consistent with the hydrostatic scaling laws (only slightly 
underdense). When t c is comparable to or longer than the loop cooling time (e.g. run 5) 
the loops begin to move along nearly cyclic sequences, spanning from hot, underdense 
conditions to cool, overdense ones. Note, however, that the range of density values (in 
dex) is significantly smaller than that covered in temperature. Moreover, the loops appear 
significantly overdense only at temperatures well below 1 MK. 

In the localized heating cases (middle panel in Fig. [3]), cadence times comparable to 
or longer than the loop cooling time (e.g. run 17) give rise to a loop evolution in general 
agreement with that found in the quasi-uniform cases with similar cadence, although the 
plots reproducing the sequences exhibit some significant differences, particularly at the 
higher temper a tures . The hydrodynamic behavior is very similar to that described by 



Spadaro et al. 



(120031 ). who considered localized transient heating on timescales comparable 
to the loop cooling time. Note that no condensation forms near the loop apex in these 
cases. When the localized heating is steady (run 6) or impulsive with t c < r coo i (e.g. 
run 9), the evolution is considerably different: as the temperature decreases, the plasma 
continues to evaporate and the density continues to rise, overshooting by a factor of ~ 5 
the equilibrium values predicted for static loops with the same coronal temperatures, until 
a plasma condensation forms near the loop apex. The fall of material on to the less heated 
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loop footpoint causes an abrupt decrease of the average density, although the loop remains 
slightly overdense and subsequently a new cycle of condensation formation and evolution 
starts. Hence these heating conditions can give rise to warm (T ~ 1-2 MK), overdense 
loops with lifetimes of some hours, such as those observed by TRACE and SOHO. 

The bottom panel in Fig. [3] shows that increasing the amount of energy supplied by 
each pulse, even without changing the other characteristics of the heating regime (e.g. 
runs 10 and 11), causes the loops to evolve toward a hot (T > 2-3 MK), quasi-static, 
slightly overdense state, where they settle after a balance between the energy supplied to 
the loop top and the radiative losses therein is achieved. Note that run 11 might apply 
to the case of the hot (~ 3 MK) loops seen by the Yohkoh Soft X-ray Telesc ope, without 



a corresponding warm counterpart (~ 1 MK) observed by TRACE (see, e.g. 



Nitta 



2000 



Klimchuk 



20061 ). If these loops are impulsively heated, then the nanoflares must occur 



frequently enough that the plasma does not have time to cool to TRACE temperatures. 



3.2. Differential Emission Measure 



Comparison of our modeling with o 



measure, DEM(T) (e.g. 



jserva tions is done using the differential emission 



Craig fc Brown 



19761 ) . which effectively describes the plasma 



distribution in temperature. 



We simulated the multi-strand loop DEM by averaging instantaneous DEMs calculated 
at n — 300 different times, randomly selected throughout the simulation. Although each 
simulation represents the evolution of a single magnetic strand, we assume that the states 
of the model at n randomly selected times can be used to describe the behavior of n 
independent strands observed at t he same time, thus giving a sin gle simulated snapshot 



of a multi-stranded loop (see, e.g. 



Patsourakos fc Klimchuk 



20051 ). The theoretical single 
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strand DEM is computed following iPeter et al.l (120061 ) 



In what follows, the obser v ed qu iet Sun (QS) and active region (AR) DEMs are 



adopt ed from 



Lanzafame et al. 



(120051 ) (SOHO-CDS observations) and 



Lanzafame et al 



(120021 ) (SERTS-89 observations), respectively. 



3.2.1. Initial conditions and hydrostatic models 

In Fig. H] we plot the DEM corresponding to: i) the initial steady equilibrium conditions 
adopted in our simulations (T apex ~ 0.75 MK, see Sect. |5J); ii) the case of a spatially uniform 
and temporally constant background heating with -Ebase = 5 x 10 -4 erg cm -3 s -1 producing 
a steady loop with an apex temperature of about 2 MK at the end of the relaxation phase 
described in Sect. |2j iii) the quasi-uniform, steady case with <E>= 4 x 10 21 erg s -1 (run 1). 
For c omparison, the ana lytical DEMs computed according to the hydrostatic equilibrium 



laws (IRosner et al 



19781 ) are also shown. The divergence of the analytical results close 
to the top of the loop comes from the DEM formula, which is written in terms of the 
inverse of a temperature gradient that vanishes at the loop maximum temperature. The 
figure shows that our quasi-uniform, steady simulation for a multi-strand reproduces quite 
closely the steady equilibrium single loop structure obtained at the end of the relaxation 
phase with a spatially uniform and temporally constant background heating of the same 
level. Apex temperatures and DEM minima at the end of the relaxation phase and in the 
quasi-uniform, stea dy simulations are essentially the same as those calculated using the 



Rosner et al. 



(119781 ) relationships. 



-19- 




Fig. 4. — Differential emission measure (DEM) for: the initial steady equilibrium conditions 
adopted in our simulations with T apex ~ 0.75 MK (thick solid line); the case of a spatially 
uniform and temporally constant background heating producing a steady loop with an apex 
temperature of ~ 2 MK at the end of the relaxation phase (dotted line); the quasi- uniform, 
steady case with < E >= 4 x 10 21 erg s _1 (run 1, thin solid line). F or comparison, the 
analytical DEMs according to the hydrostatic scaling laws (see Eq. C4 in 
are also shown (heavy- and light-dashed lines, respectively). 



Rosner et al. 



19781 ) 
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4.5 5.0 5.5 6.0 6.5 7.0 

LOG (T [K]) 



Fig. 5. — DEM from quasi- uniform impulsive simulations: run 2 (dashed line); run 3 (dash- 
dotted line); run 4 (dotted line); run 5 (dash-triple-dotted line). Comparison is done with 
QS (gray thick dashed line), AR (gray thick solid line), the quasi-uniform steady case (run 1, 
solid line), and a single pulse simulation (long-dashed line). 
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3.2.2. Quasi-uniform impulsive cases 



DEMs obtained from our quasi- uniform impulsive simulations are shown in Fig. [5j 
Since no condensation forms, the DEM is controlled essentially by the mean energy 
dissipated per unit time. The runs reported in Fig. have fixed Ep = 10 24 erg, and 
therefore the shorter t c , the higher < E >. The < E > values chosen in Fig. [5] correspond 
to cases in which the models reproduce either the observed AR DEM around logT = 5.5 
(and the QS DEM peak at logT rj 6), or the AR DEM temperature peak at logT f=a 6.3, 
plus a couple of cases with energies in between. Also shown in Fig. [5] a simulation with a 
single pulse with Ep = 10 24 erg, which reproduces very closely the conditions explored by 



Klimchuk et al. 



(120081 ). From the comparison with the observed DEMs, we see that these 
models are unable to reproduce both the TR and coronal DEM, neither for the QS nor for 
the AR. 

The quasi-uniform impulsive model with the highest mean energy dissipated per unit 
time (run 2, t c = 250 s, < E >= 4 x 10 21 erg s _1 ) is compared with the steady- uniform 
model with the same mean dissipated energy rate (run 1), confirming that there is no 
appreciable difference between the two. 

Note also that by increasing t c well above r coo i (e.g. run 5 with t r = 2000 s), we 



obtain a DEM structure very similar to that presented by iKlimchuk et al. 



(J2008|). A close 



comparison with their simulations has been made by calculating the DEM for a single pulse 
and spatially uniform heating. The res ulting DEM, also shown in Fig. [5j is almost identical, 



above logT = 5.5, to those in Fig. 2 of 



Klimchuk et al. 



(120081 ) 
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LOG (T [K]) 



Fig. 6. — DEM for impulsive localized models with different combinations of t c and Ep 
giving the same mean energy dissipated per unit time <E>= 4 x 10 21 erg s _1 : run 9 (solid 
line); run 14 (dashed line); run 16 (dash-dotted line); run 18 (dotted line). Comparison is 
made with the AR DEM (gray thick solid line). 
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3.2.3. Localized impulsive models 

In Fig. [6] we consider DEMs obtained for impulsive localized models with different 
combinations of t c and Ep giving the same mean energy rate. Choosing < E >= 
4 x 10 21 erg s _1 , we obtain a DEM peak similar to that observed in the AR. 

The mean dissipated energy determines the temperature at which the DEM is at its 
maximum (T peak ). Increasing t c up to t c R3 r coo i causes an increase in the DEM above T peak 
and a little change in the DEM slope below T pea k. Increasing further t c (t c > r coo i), causes a 
further increase in the DEM in the high temperature range (above logT fa 6.5 in the cases 
shown in Fig. [6]) and the disappearance of the coronal DEM peak. 

The differences between the cases in which condensation occurs (e.g. run 9) or does 
not occur (e.g. run 14) are too small to be appreciated in practice. Also, the localized 
steady model DEM (run 6) is almost identical to that of the localized impulsive model 
with the same <E> and t c = 250 s < r coo i- Condensation does not produce appreciable 
differences in the DEM because the bulk of the plasma maintains a temperature above 
~ 1 MK throughout the evolution (compare with Figs. [3] and [I]). The sensitivity to t c , on 
the other hand, derives from the fact that, when t c increases well above t coo i, condensation 
does not occur anymore and the bulk of the plasma has a temperature oscillating over a 
rather large range. In our simulations, such oscillation produces a smearing of the DEM 
around ~ 1 MK and eventually leads to the disappearing of the DEM peak. 

Also in this case, the simulations are unable to reproduce the whole observed DEM 
structure. By appropriately selecting t c and Ep, it is possible to reproduce the DEM in the 
high temperature range or around the minimum DEM, but it is not possible to reproduce 
the DEM in both ranges with a single (t c , Ep) pair. We conclude, therefore, that the 
assumption of multi-strand structures subject to localized impulsive heating is not sufficient 
to explain the well known discrepancies with observations, which must be due to physical 
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proces ses not included in our simulation. The recent suggestion put forward by [Judge 



( 120081 ) of cross-field diffusion of neutral atoms from cool threads extending into the corona 



may help in solving problems like this. Spicu 



below 1 MK, as suggested by 



Klimchuk et al- 



ar ab sorptio n of some of the 



fl2008h and 



De Pontieu et al. 



asma emission 



fj2009h . might 



also contribute to explain the discrepancies with the observed DEM structure. 

Nevertheless, we suggest that, despite such discrepancies, the systematic behavior of 
the simulated DEM with t c would maintain its validity even if the models are not sufficiently 
detailed, and would help in discriminating at least among very different conditions. In fact, 
despite the small sensitivity of the DEM to variations in t c , when this increases well above 
t coo i, the coronal DEM peak tends to disappear. Such changes in the DEM shape can be 
reliably verified by comparison with observations, since they are well above the expected 
uncertainties in the DEM reconstruction. The very existence of the coronal DEM peak, 
therefore, is an indication of the existence of heating pulse cadence time shorter than r coo i- 

Finally, in Fig. [7] we show the simulations with <E>= 5 x 10 20 erg s _1 (runs 7, 12, 
and 17) compared with the observed QS DEM. The dependence on t c described above for 
the AR DEM is maintained at such lower energy, which corresponds to a DEM peak close 
to the observed QS one. The DEM peak decreases at increasing t c , tending to disappear for 
t c > Tcooi- The DEM at temperatures above T pea k increases with increasing t c . 

Comparing the simulated DEM's peaks with AR and QS observations (Figs. [6] and [7]), 
the pronounced peak in the observed AR DEM may indicate a predominance of t c < r coo i 
conditions. On the contrary, the less pronounced peak in the observed QS DEM may 
indicate that t c > t coo \ conditions dominate in that regime. Although the shape of the 
observed DEM may be affected by the smoothing imposed by the regularization technique 
used to reconstruct the DEM from the observed spectral line intensities, the differences 
implied by the models are higher than the uncertainties in the DEM reconstruction due to 
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Fig. 7. — DEM for simulations with < E >= 5 x 10 20 erg s 1 (run 7, solid line; run 12, 
dashed line; run 17, dash-dotted line) compared with the QS DEM (gray thick solid line). 
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regularization smoothing and intensity ratios of spectral lines formed around and above the 
DEM peak could be useful to discriminate among the high- and low-cadence regimes. 

It is also worth noting that, despite a close fit to the observed DEM is outside 
the scope of this work and only a general comparison is made, models that reproduce 
more closely the AR DEM have a coronal electron pressure logP P <y 15.3 (P e in units 



cm 3 K) remarkably close to the value deducted by 



Lanzafame et al. 



(120021 ) using line-ratio 



diagnostics (logP c ~ 15.8). In the QS case, models that reproduce more cl osely the observed 



. DEM have logP e ~ 14.0 in the coronal part, in close agreement with the lLanzafame et al. 



((20051) estimate (logP e ~ 13.7 - 14.0). 



4. Conclusions 

We have carried out hydrodynamic simulations of multi-stranded coronal loops with 
a code capable of resolving the transition region sections and following their evolution in 
order to outline signatures of the heating regimes in the observed DEM. 

Quasi-uniform or localized impulsive heating with a cadence time ~ 1/4 of the plasma 
cooling time can produce a plasma evolution which is essentially indistinguishable from the 
corresponding steady cases. 

Plasma condensation occurs on a limited range of heating parameters. Energy 
localization is necessary to yield a catastrophic cooling phase during the loop evolution, but 
it is not always sufficient. The crucial point is the balance between the energy supplied to 
the loop top and the radiative losses therein. Thus, the variation of the energy deposition 
parameters (such as the pulse cadence, the pulse energy, or the heating damping-length) 
with respect to the global characteristics of the model (the radiative cooling time or the 
loop length) could, in some cases, prevent the occurrence of a dynamic cycle of plasma 
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condensation formation even in the presence of a localized heating. This may happen, for 
instance, by increasing the heating rate at the loop footpoints or the ratio of the heating 
damping-length to the loop length, or considering nanoflare cadence times comparable to 
or longer than the characteristic radiative cooling time. 

The DEM is found insensitive to the presence of condensation because the sequence of 
catastrophic cooling, draining toward the loop footpoint, reheating and evaporation does 
not effectively redistribute the plasma over temperature and the global distribution remains 
very close to the corresponding steady configuration. On the contrary, pulses with cadence 
longer than the plasma cooling time produce temperature oscillations in the bulk of the 
plasma which effectively smear the coronal DEM structure. The effects are observable since 
the coronal DEM peak tends to disappear when the pulse cadence is about 2 times the 
plasma cooling time. The pronounced DEM peak observed in active regions would indicate 
a predominance of conditions in which the cadence time is shorter or of the order of the 
plasma cooling time, whilst the structure of the quiet Sun DEM suggests a cadence time 
longer than the plasma cooling time. 

The warm over dens e and hot underdense loops o b served by TRACE, SO HO 



( lAschwanden et al 



1999, 



2001 



Yohkoh (jPorter fc Klimchuk 



Winebarger et al. 



2003 



Patsourakos et al. 



200J) and 



19951 ) could be explained by the dynamic evolution of the 
plasma. In particular, a localized heating producing plasma condensation cycles can give 
rise to warm (T ~ 2 MK), overdense loops with lifetimes of some hours, such as those 
observed by TRACE and SOHO. Moreover, increasing the amount of energy supplied by 
each pulse at the loop footpoints prevents the plasma condensation formation and causes 
the loops to settle in a hot (T ~ 2-3 MK), quasi-static, slightly overdense state, which 
can reproduce the case of hot loops se en by Yohko h SXT, without a corresponding warm 



counterpart observed by TRACE (e.g., 



Nitta 



2000) 
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Nevertheless, our simulations are unable to reproduce both the transition region and 
the coronal DEM structure with a unique set of pa ramete r s, wh ich sug gest that some 



additi onal physical processes , like t hat proposed by 



fl2008h and 



De Pontieu et al. 



Judge! fj2008h or by 



Klimchuk et al 



(120091 ). must be taking place in the transition region. 
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